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ABSTRACT 

New infrared absorption measurements of oxygen isotope ratios in CO gas from individual young 
stellar objects confirm that the solar system is anomalously high in its [^^0]/[^^0] ratio compared 
with extra-solar oxygen in the Galaxy. We show that this difference in oxygen isotope ratios is best 
explained by ~ 1% enrichment of the proto-solar molecular cloud by ejecta from type II supernovae 
from a cluster having of order a few hundred stars that predated the Sun by at least 10 to 20 Myr. 
The likely source of exogenous oxygen was the explosion of one or more B stars during a process of 
propagating star formation. 

Subject headings: Galaxy: abundances — infrared: stars — radio lines: ism — Sun: abundances 



1. INTRODUCTION 

Solar system oxygen isotope ratios are peculiar with 
respect to Galactic values. All solar system materials 
have [i*0]/["0] ratios of 5.2 ± 0.2 ([^^O] refers to the 
abundance of ^*0 by number) (see Appendix). Values 
from measurements of radio emission from isotopologues 
of CO, OH, H2CO and HCO"*' in molecular clouds across 
the Galaxy, on th e other han d, are between 3.5 ib 0.3 
(pannier 1989; P enziasj [l98T| and 4.1 ± 0.1 ( [Wouter-^ 
pflot et al.||2008 |). I'igure [l] shows radio emission data 
tor CO oxygen isotopologues in molecular clouds span- 
ning a large range of distances from the Galactic center. 
Both the older data ( Penzias|1981 ) and the revised data 
(Wouterloot et al. 2008p are included. The latter differ 
from the former by inclusion of more rotational transi- 
tions and a model for optical depth effects. Also shown 
is the [*^0]/[^^0] for the solar system. The radio data 
are shown together as probability density contours that 
are essentially a smoothed two-dimensional histogram. 
The method for contouring the data is described in §4.2. 
It is clear from Figure [T] that the solar ratio of the rare 
oxygen isotopes is inconsistent with the vast majority of 
Galactic values. The newer data also suggest a hint of a 
trend towards greater ^^O relative to ^®0 in the Galactic 
center and an excess of ^^O relative to "'^''O in the outer- 
most Galaxy. We note that in view of radial rates of 
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Fig. 1.— Plot of [Ci*0]/[C"0] vs. distance from the Galactic 
center (Rqc;) for previously published molecular cloud data | |Pen-| 
zias 1981 Wouterloot et al. 2008) (contours), the solar system, and 
the new survey ol young stellar objects (data points with 2cr error 
bars). The intensity of grey scale from light grey to black shows 
the probability density of previous molecular cloud [C^*0]/[C^'^0] 
measurements. 



stellar migration in the Galaxy, the Sun could not have 
formed near the Galactic center nor could it have formed 
beyond ^1 kpc from the Galactic center (e.g 
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Explanations for this disparity in solar and Galactic 
oxygen isotope ratios have included: 1) systematic er- 
rors between radio emission data used to measure molec- 
ular cloud f^®0]/[^''0] and m ass spectrometry meth ods 



used to measure solar values (Prantzos et al. 1996); 2 



a burst of high-mass star formation produced by merger 
of another galaxy with t he Milky Way approximately 5 
to 6 Gyr before present ( Clayton||2003 1 ; and 3) Galactic 
chemi cal evolution o f ["'"^OJ/["^ ''OJ over the past few billion 
years (|Nittler||2009|). In this paper we suggest that the 



peculiar [^*0]/[^'0] of the solar system is a consequence 
of local enrichment of the solar birth environment by ex- 
plosion of nearby B stars. The work is motivated by our 
recent infrared abso rption data for young stellar objects 
(jSmith et al. 2009) showing that systematic errors are 
not likely to be the cause of the disparity. We show that 
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the oxygen isotopic composition of the exogenous type 
II supernova ejecta required to explain the solar ratio of 
rare oxygen isotopes places constraints on the maximum 
masses of the supernova progenitors. The limit on pro- 
genitor masses in turn suggests a scenario of enrichment 
by propagating star formation leading up to the forma- 
tion of the Sun. This scenario is distinct from previ- 
ous suggestions of supernova enrichment in which it had 
been assumed that t he supernova progenitors were coeval 



with the Sun (e.g., Schramm & Olive 1982). It is con 



sistent, however, wixh a proposed origin tor short-lived 
radionuclides in the early solar system in which several 
supernovae from a previous generation of star formation 
enrich nearby molecular cloud ma terial prior to collaps e 
to form a new generation of stars ( [Gounelle et al.|[2009 ) . 

The isotopic compositions of oxygen are reported here 
in plots of [i^O]/[i^O] against [180]/[160], referred to 
commonly as "three-isotope" plots, throughout this pa- 
per. Because the differences in isotope ratios we are 
concerned with are large, we will use a variant on the 
delta notation commonly used in the cosmochemistry 
literature to report isotope ratios. In this notation, 
(5170' = 103 ln([i701/[i60],/p70]/ roTTo„ ) where ISM 



refers to the local ISM as defined in Wilson ( 1999 ) (an ar- 
bitrary but convenient reference with~p^O]/p^O] = 557 
and [i'^0]/["0] = 2005.2) and i refers to an object or 
mass interval of interest. Positive values for J^^O' sig- 
nify enrichment in ^^O relative to ^^O with respect to the 
local ISM while negative values for 6^^0' signify a rela- 
tive depletion. Values for S^^O' are defined analogously 
(the prime in each instance signifies this logarithmic form 
of the delta notation) . 

2. GALACTIC CHEMICAL EVOLUTION OF i^O/i^O 

Both i''0 and ^^O are secondary nuclides, produce d 
by H and He burning, respectively (Meyer et al. 2008). 
One expects their solar ratio to be that ot the bulk of the 
Galaxy at the time of solar system formation unless the 
abundances of the oxygen isotopes were affected by some 
local perturbation. Quantitative models for the increase 
in the abundances of the oxygen isotopes with time as 
part of the overall Galactic chemical evolution (GCE) 
show that we should expect that even as [^^0]/[^^0] and 
[^*0]/[^®0] rose linearly with time, the ratio of the two 
secondary nuclides, [^^0]/[^''0], shou ld have been nearly 
constant after the first bill ion years (Woosley & Weaver 
1995[ [Prantzos et al.|[T996l ). A constant ['*OJ/[^'OJ over 
time would mean that the anomalous value for the solar 
system compared with the present-day Galaxy cannot be 
attributed to GCE over the past 4.6 billion years (the age 
of the solar system). 

A trend of decreasing [^^0]/[^^0] and [i^O]/[i^O] ra- 
tios at nearly constant [^*0]/["'^^0] with increasing galac- 
tocentric radius, R gc , was well e stabhshed by early ra- 
dio emission data (Wilson 1999|, and the most recent 
data confirm this trend (l^'igures [2| and [3| . These trends 
of varying oxygen isotope ratios with Rgc ai'e consistent 
with our expectations from GCE, where Rgc is used as 
a proxy for time (a well-established means of investigat- 
ing chemical evolution in the Galaxy). For this reason 
we conclude that it is unlikely that a shift in interstellar 
[^*0]/[^'^0] from 5.2 to 4.1 occurred over the last half to 
third of the lifetime of the Milky Way at plausible loca- 
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Fig. 2.— Plot of [i8o]/[l6o] vs. [i'^0]/[i'*0] for molecular clouds 
across the Galaxy in delta notation. Molecular cloud data (squares) 
represent I'^CO, C^*0 and C^^O abundances from Wouterloot et 
al. (2008) combined with best-fit I^CO/^^CO vs. Rqc data from 
Milam et al. (2005). The effect of using the [i^q/fi^C] vs. Rqc 
data for H2CO and CN rather than CO presented by Milam et 
al. (2005) is to move the molecular cloud points up and down the 
slope-1 line by approximately 200 per mil, but the slope-1 line is 
preserved. A best-fit line to the molecular cloud data with equal 
weighting yields a slope of 1.11 it 0.04 where the la uncertainty 
applies assuming ±100 for absissca and ordinate for each datum 
(consistent with a reduced of unity). Also shown are measure- 
ments for the solar system compiled from the literature (circles) 
and 95% confidence error ellipses for the young stellar objects from 
Smith et al. (2009) also shown in Figure [l] 




Fig. 3.— Plot of [Ci*0]/[C"0] vs. distance from the Galactic 
center (Rgc) expressed in delta notation relative to the local ISM. 
Data represent ^^CO and C^*0 abundances from Wouterloot et 
al. (2008) combined with best-fit ^^CO/^^CO data from Milam et 
al. (2005). 



tions for the formation of the Sun 
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( 2009 ) offer an alternative view of Galac- 
tic chemical evolution of oxygen isotopes. In that work 
the authors present results of a two-box model calcula- 
tion in which they posit that the contribution of ^^O from 
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AGB stars results in a substantial decline in [^^0]/[^^0] 
with time in the Galaxy. Details of the calculations (e.g., 
equations representing the box model) are not presented, 
but the authors note that their model predicts differences 
between IS M and s tar forming regions that arc not ob- 
served (e.g. Smith et al.^2009J . They also note that their 
model is unable to reproduce the oxygen abundances of 
the ISM 4.6 Gyr bef ore present, presuma bly meaning 
the solar value. The [Gaidos et aL (2009) model pro- 
vides no explanation tor the trends m oxygen isotope 
ratios shown in Figures [2] and |3] Their GCE evolution 
model has a slope of nearly 3 on an oxygen three-isotope 
plot rather than the slope of 1 defined by the molecu- 
lar cloud CO isotopogue data in Figure [2j The calcula- 
tions were meant to explain the difference between local 
ISM and solar ['^^0]/[^'^0], but other than this discrep- 
ancy (the subject of this paper) we are aware of no data 
representing the Galaxy that define a slope of ~ 3 in 
three-isotope space. In support of their model the au- 
thors suggest that there is no evidence for the slope- 1 
line in oxygen three-isotope space. We find this asser- 
tion difficult to defend in view of the CO data shown in 
Figure [2j Their assertion is made wi th reference to seven 
measurements of P^OH]/[^*OH] by Polehampton et al. 



( 2005| ) that show a mean [^^OR]/[^^0}i] of 424 ± 70 la 
(corresponding to a S^^O' value of -t-285 ± 168) and no 
apparent trend with Galactocentric radius. This result 
is at odds with the substantial gradient shown in Fig- 
ure |3] based on the most recent measurements of molec- 
ular cloud CO isotopologues. Reasons for there being 
a twenty fold increase in [Ci60]/[Ci*0] wit h Rqc but 
nearly constant [^^OH]/[^^OII] are unclear. Polehamp- 
ton et al. (2005) point to low number statistics tor their 
OH measurements and lack of constraints on some of the 
Galactocentric distances used in their work. In any case, 
because CO is the dominant gas-phase reservoir of oxy- 
gen in molecular clouds (e.g., Wakelam et al.pOlO ), we 
expect this molecule to be, arguably, the most robust 
tracer for oxygen isotope ratios across the Galaxy. 

3. OXYGEN ISOTOPE RATIO DATA FOR YOUNG 
STELLAR OBJECTS 

The possibilities for systematic errors in the radio data 
for molecular clouds, or biases due to disparate scales 
of observation for clouds and the solar system, are ad- 
dressed by new data for young stellar objects (YSOs). 
We use results from an infrared (IR) spec troscopic sur- 
vey of mole cules in young stellar objects (Pontoppidan 
et al. 20081, part of which aim s to measure oxygen iso- 
tope ratios with h igh precision ( [Pontopiddan et al.|2008 
Smith et al.|2009|). These new data comprise IK absorp- 



tion spectra tor CO rather than radio emission and rep- 
resent scales of observation commensurate with that of 
the solar system (on the order of hundreds of astronom- 
ical units). By embodying an entirely different scale of 
observation and an independent method of measurement 
these new data should circumvent sampling bias and/or 
systematic errors that might be present in the radio emis- 
sion results. The purpose in this context is to establish 
the magnitude and scale of oxygen isotope variability in 
young stars for comparison with the solar system. Our 
survey now includes column densities for CO isotopo- 
logues of gas surrounding three young stellar objects, in- 
cluding Reipurth 50 (RE 50), VV Corona Australis (VV 



CrA), and IRAS 19110-^1045. High-resolution 4.7 ^lm 
fundamental and 2.3 /xm overtone ro-vibrational absorp- 
tion bands for CO were obtained for RE 50 and VV CrA 
using Cryogenic High-Resolution Infrared Echelle Spec- 
trograph (CRIRES) on the Very Large Telescope (VLT) 
at ESO's Paranal Observatory. Lower spectral resolu- 
tion data from Near Infrared Echelle Spectrograph (NIR- 
SPEC) on the Keck II telescope were also analyzed for 
a more massive embedded young stellar object, IRAS 
19110+1045. Details of the data collection and process- 
ing for RE 50 and VV CrA are provided by [Smith et al.| 
(j2009) as part of their study of the potential for pho- 
tochemical effects on the isotopic composition of CO in 
YSOs . Thos e for IRAS 19110-^1045 are given by [Smith 
eFdl ( [2007| ). RE 50 is an embedded YSO in the Orion 
star-tormmg cloud ~ 470 pc from the Sun. It is an 
FU Ori type object in stage I. VV CrA is a stage II 
T-Tauri disk 130 pc from the Sun. The third object, 
IRAS 19110+1045, is a more massive embedded YSO at 
a Galactocentric radius, Rgc, of 6 kpc. 

The YSO data are compared with the molecular cloud 
radio emission data in Figures fl] and [2] These results 
are consistent with the molecular cloud radio emission 
data in showing that typical Galactic [^^0]/[^^0] is near 
4. The [Ci**0]/[C"0] values are 4.1 ± 0.4 for VV CrA, 
4.4 + 0.2 for RE 50, and 4.0+1.7 for IRAS 19110+1045. 
We conclude from these measurements that the solar sys- 
tem is indeed unusual, and that the difference between 
solar [180]/[170] and typical Galactic [i*^0]/[i^O] can- 
not be attributed solely to heterogeneity in oxygen iso- 
tope ratios on the scale of individual stars. A systematic 
error in the molecular cloud data is also now excluded. 

4. SUPERNOVA ENRICHMENT OF THE PROTOSOLAR 
MOLECULAR CLOUD 

4.1. The Oxygen Isotopic Compositions of Supernovae 

We find that the most likely explanation for the anoma- 
lous [^®0]/[^''0] of the solar system is that the Sun and 
planets formed from molecular cloud material enriched in 
exogenous ^^0-rich oxygen ejected from stars in a nearby 
cluster that ended their lives as type II supernovae (SNe 
II) . Measurements of the relative abundances of all three 
stable oxygen isotopes place critical constraints on mix- 
ing between normal Galactic oxygen and this ^®0-rich 
ejecta. We explore this proposition, and its consequences 
for the origin of the solar system, below. 

With three isotopes of oxygen we can compare the 
isotopic composition of the ISM at the time the Sun 
was formed with different supernova products in order 
to identify plausible sources of exogenous oxygen. The 
oxygen isotopic composition of the local ISM 4.6 Gyr be- 
fore present can be estimated from ages of between 13.6 
Gyr and 10 Gyr for the Galaxy and the 4.6 Gyr age of 
the solar system. As described above, both [^*0]/[^^0] 
and [^^0]/["'^^0] are ratios of secondary to primary nu- 
clides that have risen linearly, to first order, with time in 
the Galaxy at nearly constant [^^0]/['^'^0]. The rise m 
the abundances of ^^0 and ^^O relative to that of ^^O 
linearly with time leads to the relationship Aage/age ~ 
(A[i80]/p60])/([i«0]/p'^0]). We therefore expect both 
[i80]/[i60] and ["0]/[i60] ratios in the ISM to have 
risen by between 35% and 46% in the past 4.6 billion 
years. More complicated models for oxygen GCE are 
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consistent with this estimate ( Prantzos et al.|1996 1 . The 
^ 350 to 460%o increase in both ratios can be subtracted 
from the composition of present-day ISM to obtain an 
estimate of the oxygen isotopic composition of the ISM 
at the time of the formation of the Sun. The precision 
of this estimate is limited by an uncertainty of at least 
-20% (±200 %o) in [i80]/[i60] and [^^0]/[^<^0] for the 
present-day ISM. With this uncertainty, the present-day 
nomi nal \^^0]/\^^0 ] and [^''0]/[^^0] values of the local 
ISM (Wilson 1999) are reasonable upper limits for the 
values tor the ilSM 4.6 Gyr ago. 

Previous work has emphasized that [^^0]/[^^0] of oxy- 
gen liberated by the explosion of a SN II varie s system- 
atically with the mass of the progenitor star ( Gounelle 
&: Meibom 2007). One can illustrate the inputs of in- 
dividual stars to the composition of stellar ejecta from 
star clusters using mass fractions over small intervals of 
the initial mass function. Starting with the initial mass 
function, yielding the number of stars N of mass m 



am 



(1) 



integration over some mass interval m; to m„ (lower to 
upper) yields for number of stars 



N„ 



and for mass 



— f3 m "dm 
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The fraction of mass contained in stars of masses m; to 
TO,, is then 



X„ 
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(4) 



From these mass fractions of the total stellar system we 
obtain mass fractions of oxygen ejected from a genera- 
tion of star formation by fitting existing models for mass 
loss as functions of progenitor mass. For stars > 8Mq 
we used the SNe II yields (by yields we mean production 
and not net yield that results from con sidering the nu- 
clides remaining in the stellar core) from Rauscher et al. 
([2002) (RHIIW02, including supplemental tables avail- 
able on line) taking into account also t he calculations of 



Woosl ey & Weaver 
|2007| (WHOTH^ 



( |1995| (WW95) and |Woosley fc Heger 
d and mass cut tables kindly provided 
by A. Heger). These studies exhibit large differences in 
^^O production. A reduction in ^^O yield from WW95 to 
RHHW02 ref lects the revision to th e destruction rate for 
that nuclide ( Blackmon et al.|1995 1 . A further reduction 
in the^^O yield from RHHW02 to WHO? resuhed from 
lower initial CNO abundances of the progenitor stars. 
The mass of oxygen ejected from the larger progenitors 
(> 3OA/0) depends on the choice of the intensity of the 
SN "piston" employed by RHHW02. For the calcula- 
tions presented here, we made use of the higher-energy 
piston models. The resulting yields are similar to those 
of WW95. The oxygen yields as a function of progen- 
itor mass are shown in Figure [4] along with the fits to 



the calculations used in the present calculations. The fit 
to the RHIIW02 calculations for mass of oxygen ejected 
relative to mass of progenitor, Mq/M*, is 



Mo 



4 l\j2 



: 1.0250 X IQ-^'M^ + 5.4290 x lO^^Af, 
1.6617 X 10"^ Af* - 1.1800 x 10"^ 



(5) 

where masses are in solar units. For comparison with the 
SNe II calculations, and for calculating mass fractions of 
oxygen relative to all of the oxygen injected into the ISM 
by a stellar cluster, we include the oxygen produced by 
mass loss from AGB stars for ikL < SMrr) as calculated 



by [Karakas fc Lattanzio] ( [2007| ) (KL07) (Figure [I]). The 
equation tor the tit to the KL07 calculations for mass of 
oxygen released relative to mass of progenitor is 



Mo 



5.5687 X IQ-Hl^ + 8.4985 x 10"*A1 



4 /l,f2 



3.9956 X lO^^Af* + 1.3560 x 10 
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(6) 

From equations (4) through (6) the fraction of oxygen 
attributable to a progenitor mass Af* in the total ejecta 
from a generation of stars, ^o,Am, represented by a spec- 
ified initial mass function is 



Mo 
M, 



Am 



(7) 



where Am represents the mass interval to; to to„, the 
summation is over all mass intervals, and Af* is the mean 
mass for that interval. We are concerned here with both 
the total oxygen ejected and the isotopic composition of 
that oxygen. Oxygen isotope abundances for the ejecta, 
expressed as the mass fraction of oxygen for nuclide i, 
X*0 = M.o/{Mi6o+Mi7o + Miso), are shown in Figure 
[5I The fit to the RHHW02 calculations used here for ^^O 
is 



0.0048 



(1 + exp(-(Af* - 18.3973)/(-2.1545))) ' 



(8) 



For ^''O we fit the RHHW02 calculations and then scaled 
the results by a factor of 0.5 to bring them into line with 
the more recent results from WII07. The fit prior to 
scaling (Figure [S]) is 



: 3.1995 X 10^'' - 1.9676 x IQ-^Af* 
-3.043 X IQ-Hl'i. 



(9) 



Using the calculations described above we can plot 
the isotopic compositions and mass fractions of oxygen 
ejected from stars comprising a stellar cluster in three- 
isotope space (Figure |6]). The result shows that SNe 
II produce a wide range of compositions from the high 
[i80]/[iS0]- low [i^O]/[i^O] ejected by smaller progeni- 
tors with masses less than 30AfQ to the low [^^0]/[^^0]- 
high [^''0]/[^^0] ejected by the more massive progeni- 
tors. It is clear that exogenous oxygen from low-mass 
SNe II (< 2OAf0) could explain the anomalously high 
[^*0]/[^''0] of the Sun compared with the more normal 
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Fig. 4. — Predicted masses of oxygen ejected by type II super- 
novae as a function of progenitor mass (Mq) as given by RHHW02 
and WW95. Also showrn are the predicted masses of oxygen re- 
leased by AGB stars as given by KL07. Fits to these calculations 
used in the present study are shown as solid lines (Equations (5) 
and (6)). Two sets of calculations by RHHW02 are shown based 
on low and high-energy piston models. 
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Fig. 5. — Mass fractions of oxygen isotopes 



MiQl {Mwq + Mitq + Mi8q), as functions of progenitor mass 
(MQ)as predicted by RHHW02. Fits used in the present study are 
shown by the solid curves. 

compositions of the Galaxy. The implication is that the 
one or more stars that exploded and enriched the pro- 
tosolar cloud with [^^0]/[^^0]-rich oxygen were B stars. 
The more massive O stars, on the other hand, eject oxy- 
gen with [^^0]/[^^0] too low to allow for mixing with 
ancient ISM to produce the solar oxygen isotope ratios 
while AGB stars, the sources of pre-solar grains in mete- 
orites, produce oxygen too low in [^^0]/[^^0] to explain 
the solar values. 

Supernova yield calculations are inherently uncertain 
and this uncertainty is difficult to quantify. In order t_o 
illustrate the level of uncertainty, we show in Figure 
examples of ejecta oxygen isotope ratios as functions 
SNe II progenitor masses from a variety of recent stud- 
ies. In all cases one sees that the lower-mass progen- 
itors in all of the calculations do indeed tend to pro- 
duce higher [180]/[160] and [i^0]/[i*50] ejecta than the 
higher mass-progenitors. However the exact locations of 
the yields relative to local ISM in three-isotope space 
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Fig. 6. — Three-isotope plot showing the isotopic compositions 
of oxygen ejected by type II supernovae (M* > 8M©) and AGB 
stars (Af, < SMq) obtained by smoothing supernova yields from 
RHHW02 with [^^O] production scaled to match that of WH07. 
Numbers adjacent circles indicate progentior masses (Mq). Each 
circle represents a progenitor mass bin of IMq. The diameters 
of the circles represent the fractions of the total oxygen released 
by all stars comprising the complete IMF that are attributable to 
the individual mass bins (Equation (7)). For this calculation a 
"generic" IMF was used where ai = 1.25 for M* < IMq, Q2 = 2.1 
for IMq < M, < 2M0, 03 = 2.3 for 2Mq < M, < IOMq, and 
a4 = 2.3 for IOMq < M*. The result is not altered significantly 
by employing a different IMF. 

vary considerably from model to model. Of the calcu 



lations considered here, those of Nomoto et al. (20061 
(NTUKM06) exhibit the greatest disparity relative to 
the others, with substantially greater abundances of both 
secondary oxygen nuclides compared with the results of 
WW95, RHHW02, Limongi fc Chiefh (2003) (LC03), and 
WH07. Our smoothed yield vs. mass tunction is consis- 
tent with the most recent of these models (WII07) by 
design (Figure [7| . The veracity of what follows depends 
on the extent to which current calculations faithfully rep- 
resent supernova oxygen isotope yields, but in all cases 
the inverse relationship between [i'*0]/[^^0] and SN II 
progenitor mass is robust. 

We note that there is a discrepancy in the produc- 
tion of 1^0 between WH07 and RHHW02. The yield 
given by WHO? is substantially greater than that given 
by RIIIIW02 for stellar masses greater than 25 Mq (Fig- 
ure n\. However, we have performed the calculations 
described below using fits to both sets of oxygen isotope 
yields and find no appreciable difference in the results. 
This is in part because we are concerned with supernova 
progenitor masses considerably less than 25 Mq. 

For the purposes of these calculations, the precise form 
of the initial mass function (IMF) is unimportant. This is 
because the main discrepancies between different IMFs 
are for the very low masses and the very high masses, 
while, as we show below, we are primarily concerned with 
the mass range from 8 to 3OM0. 

Given the apparent requirement that exogenous oxy- 
gen that affected the isotopic composition of the solar 
system came from B stars but not O stars, the ques- 
tion then arises as to how oxygen from exploding B stars 
would be "selected" by the solar precursor in greater 
abundance than oxygen from other sources. The answer 
lies in the stochastic nature of star formation in general. 
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Fig. 7. — Oxygen isotope ratios of SNe II ejecta as functions 
of progenitor mass (Mq) as predicted by LC03, NTUKM06, and 
WH07 for initial compositions of solar metallicity. The fit used in 
the present study (Figure rol is shown for comparison. Also shown 
are data representing the solar system and the position of the local 
ISM as adopted in this paper. Numbers next to each datum are 
the progenitor masses {Mq). 

as described below. 

4.2. Oxygen Isotope Ratios of Ejecta as a Function of 
Cluster Size 

The relative age of the cluster proximal to the proto- 
solar molecular cloud is limited by the constraint on the 
average supernova progenitor star mass. B stars having 
masses < 20Mq and > 8Mq (the minimum mass to pro- 
duce a type II SN) re quire 10 to 30 Myrs to evolve prior t o 
explosion as SNe II ( |Prantzos|2008l [Schaller et aL|1992] . 
The stars that were the source of the high [^'^0]/[^'0] 
oxygen therefore belonged to a generation of star forma- 
tion that predated the solar system by at least 10 Myr. 
We emphasize this conclusion. The oxygen isotope data 
indicate an episode of enrichment that involves a genera- 
tion of star formation that predated that which produced 
the Sun. This conclusion is in contrast to, and inconsis- 
tent with, earlier m odels for supernova enrich ment of the 



solar system (e.g., Schramm & Olive 1982) in which it 



is assumed that the exploding supernova and the Sun 
formed together as siblings in the same star cluster. 

Enrichment of a region of star formation by explosions 
of B stars from an earlier generation of star formation 
is consistent with protracted star formation in molecular 
cloud complexes in the G alaxy spanning 10 to 20 Myrs 
(e.g.,'Go unelle et al.||2009[ ). With this scenario of propa- 
gating star formation m mind, we used a statistical anal- 
ysis to examine the likelihood for [^*0]/[^^0] enrichment 
of the proto-solar molecular cloud by oxygen ejected from 
B stars that evolve to become SNe II . We adopted 20 
Myr as a conservative upper limit for the time interval 
over which B stars reside in a cloud complex after for- 
mation, corresponding to a minimum stellar mass of 11 
Mq. Current evidence indicates that young clusters of 
moderate size that produce one or more B stars disrupt 
their par ental molecular clouds on time scales o f ~ 3 to 
10 Myr ( [Leisawitz et al.||1989] |Elmegreen||2007[ ), a time 
span less than the 10 Myr required tor the most massive, 
and therefore the most short-lived, B stars {18Mq, ~ 



BO) to explode as type II supernovae. However, molec- 
ular clouds are usually present in cloud complexes ex- 
tending over hundreds of parsecs, and these complexes 
are sites of protracted epis odes of star formation last- 
ing ten s of millions of years ( Elmegre en|2007 Hartmann 



et al. 20011. With velocity dispersions ot ~ lU km/s 



betwe en stellar subgroups within a complex (de Bruijne 
19991, the more massive B stars (BO to Blj can drift 
only 100 to 200 pc during their lifetimes. Such a star 
cannot leave the vicinity of the giant molecular cloud 
(GMC) complex before exploding as a SN II; BO to Bl 
stars formed in one cloud are likely to be in close prox- 
imity to other clouds within the same extended complex 
when they become supernovae. Indeed, such supernovae 
are often invoked as triggers for succe ssive generations 



of star formation (Briceno et al. 2007L Evidence that 



star clusters can encounter multiple clouds is provided 
by the presence of multiple generations o f stars within 
some well studied clusters (Mackey 2009). The net re- 
sult is that complexes hundreds of pc m size persist for 
up to 50 Myr, GMCs within the complexes survive for 
10 to 20 Myr, and G MC cores t hat produce individual 
clusters last ~ 3 Myr ( jElmegreen, 2007 ). 
For the s tatistical analysis we used the m ass generation 



functi on of Kroupa et al. 



^ ^ 1993[ ) modified by Brasser et al 

( 2006 1 to obtain 300 random realizations of star clusters 
of various sizes ranging from tens to tens of thousands of 
members. With this function the mass of the jth star of 
the cluster is obtained from the expression 



Mj/M© =0.01 + (0.19x'-''^+0.05x''-'')/(l-x)"-^« (10) 

where 2: is a uniformly distributed random number be- 
tween and 1. Sampling of the IMF, as simulated with 
equation ( 10 ), results in systematic relationships between 



occurrences of SNe II, the maximum size of SNe progen- 
itors, and cluster size. Figure [8] shows that the aver- 
age fractional number of stars that produce type II SNe 
within 20 Myr of the birth of the cluster is 1.4 x 10"'^ 
for all cluster sizes. The fact that all of the simulated 
clusters produce the same average fractional number of 
supernovae illustrates that clusters of all sizes do indeed 
represent random samplings of the IMF. However, the 
discrete, stochastic sampling produces important differ- 
ences in populations as a function of cluster size. Clusters 
composed of relatively few numbers of stars can some- 
times produce no SNe II at all simply because the pro- 
genitor stars (those with > 8Mq) are relatively rare. 
This is illustrated in Figure [9] where it can be seen that, 
with the 20 Myr time constraint where we consider only 
stars with masses greater than 11 Mq, only about 15 % 
of clusters compo sed of 100 star s, a p ractical minimum 
size for clusters ( Lada fc Lada||2003 ), produce SNe II 
whereas 80% of clusters ot 1,000 stars produce SNe II 
in these simulations. The maximum stellar mass, and 
therefore most massive supernova progenitor, also varies 
with cluster size as a fundamental consequence of dis- 
crete sampling of the IMF in which the more massive 
the star the more rare its occurrence. As cluster size 
increases, the maximum in the frequency of most mas- 
sive members of the cluster shifts to higher mass (Figure 
10). In other words, smaller clusters tend to produce 
smaller supernovae and larger clusters tend to produce 
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300 draws, 20 Myr 



N(SNts) = 1.4347x10 N. - 9,509x10 




Fig. 8. — Plot of mean number of type II SNc (ordinate) occur- 
ring within 20 Myr of birtli of star cluster versus number of stars 
comprising tiie cluster (abscissa) pro duce d by 300 random draws 
for each mass interval using equation \10\ . The 20 Myr constraint 
corresponds to a minimum stellar mass or 1 1 Mq . Error bars rep- 
resent the 1(T distribution for the random draws at each cluster 
size. All cluster sizes produce the same fractional number of stars 
that become SNe II (1.4 X 10~^), showing that all of the clusters 
obey random sampling of the same initial mass function. 




Fig. 9. — Plot of fraction of clusters composed of Nt stars that 
produce SNe II within 20 Myr of formation of the cluster as pro- 
duced using equation \10\ . A fit to these simulation results pro- 
vides an expression for the the fraction of clusters that produce 
supernovae II, /, as a function of number of stars comprising the 
cluster: / = (1 - exp{-1.5A X lO-^Af,)). 

larger supernovae. 

We emphasize two important points in this regard. 
First, the histograms in Figm^e [lO] do not depict the 
frequency distributions of all masses, only that of the 
most massive member of each cluster. Second, that 
the propensity for larger clusters to more reliably pro- 
duce the more massive stars is a natural consequence of 
stochastic sampling of any initial mass function in which 
frequency varies inversely with stellar mass. Although 
debate surrounds wheth er sampling of the IM F by star 
clusters is truly random ( Weidner et al.|20T0 | , the posi- 
tive correlation between the mode in maximal star mass 
and cluster size is evident nonetheless (e.g., 
eTaLjlMol figure 4) 



Weidner 



We combined supernova oxygen isotope yields with the 
statistical analysis described above in order to examine 
the relationship between star cluster size and the oxy- 
gen isotopic composition of the oxygen ejected by SNe 
II from the cluster. The results are depicted in oxygen 
three-isotope space by contouring the relative probabil- 
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Fig. 10. — Histogram showing the frequency of occurrences of the 
maximal stellar mass produced in each of IC* random instances of 
clusters composed of 500 stars (solid) and 2000 stars (lined). Note 
that the mass of the largest star tends to be greater for the larger 
clusters. In particular, the peak in the mass of the largest star shifts 
from 8 Mq for A^, = 500 to approximately 20 Mq for Af, = 2000. 

ity of occurrence of supernova oxygen ejecta of a given 
isotopic composition. The isotope ratios of the ejecta 
are integrated from time to 20 Myr after the (instan- 
taneous) formation of the cluster. Assignment of prob- 
abilities is an expediency that amounts to a Gaussian 
smearing of each model datum. In this way clusters of 
adjacent points receive greater weight (greater probabil- 
ity) than individual points. For this purpose each of 300 
time-integrated oxygen isotopic compositions of SNe II 
ejecta is smeared by a Gaussian distribution: 



Pl70 = 

and 

Piso = 



: CXp 



1 /<5^^0'-<5^^0'Modcl 



(11) 



: CXp 



1 /<5i«o'-<5i8o'Modcr 



(12) 



where (J^^OVodci and are the model super- 

nova ejecta isotope ratios in S' notation and a is taken to 
be slightly larger than the grid spacing. The normalized 
joint probability for a given grid square in three-isotope 
space is then 



dR 



1 



grid 



(13) 



where dS^'^O' = dS^^O' = a and n is the number of points 
(300 in this case). For the results shown here we set a 
equal to a convenient but arbitrary level of smearing of 
40%c compared with a grid spacing of 33%c. Altering 
a and the grid spacing has no substantive effect on the 
results. 

Probability density contours for the isotopic composi- 
tions of oxygen ejected over a 20 Myr period from clusters 
of 500 and 5,000 stars are shown in Figure [TT] The con- 
tours quantify the inverse relationship between cluster 
size and [^^0]/[^^0] of oxygen ejected by SNe II. Oxy- 
gen ejected en masse from a cluster of 5,000 stars (grey 
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5"0' = 10'ln( ('"0/"0y("'0/"0) |g„ ) 



Fig. 11.— The ISM {[^^ O] / [^'^ O] = 4.1), Sun (circle/dot), and 
example results from the statistical analysis (contours) shown in 
oxygen three-isotope space. Black contours show the probability 
density for 300 realizations of the oxygen isotopic compositions 
of ejecta from clusters of 500 stars emitted over a 20 Myr time 
interval. Grey contours show the probability density for ejecta from 
clusters of 5,000 stars. The contour interval is 0.005 to a maximum 
of 0.04. The two mixing curves show the range of possible ways 
of mixing SNe II ejecta with ISM oxygen 4.6 Gyr before present 
having a typical Galactic [^*0]/[^'^0] of 4.1 to form the Sun and 
solar system with [l*0]/[i^O] of 5.2. 




8"0' = 10^ ln(("0/'^0)./('WO)|g„) 

Fig. 12. — Three-isotope ratio plot showing the 95% confi- 
dence ellipses for YSOs VV CrA and RE 50 and the YSO IRAS 
19110-1-1045 datum. Errors in the IRAS 19110-1-1045 data are 
dashed because they are poorly characterized in [^*0]/[^^0] vs. 
[^'^Ol/I^^O] space as a consequence of using a curve-of-growth 
analysis of limited precision to obtain the C^^O column density. 
Dots on the ISM-SNe II mixing curves represent 10% increments 
of addition of SNe II oxygen to molecular cloud oxygen with 
f-'-^Oj/f -'-^O ] = 4.1. The two curves are the same as those shown in 
Figure [TT] The steeper mixing curve corresponds to mixing with 
ejecta from an average SNe II progenitor mass of IYMq while the 
shallower curve shows mixing with debris from a SNe II progenitor 
mass of 12 Mq. 



contours) is more likely to extend to lower [i^O]/[i^O] 
than is oxygen ejected from a cluster composed of 500 
stars (black contours). 

We also show in Figure [TT] calculated mixing curves 
between estimates for the ISM 4.6 Gyr before present 
and the most likely compositions produced by a clus- 
ter of ca. 500 stars that would produce the composi- 
tion of the Sun. One curve is based on mixing with the 
most probable composition for a 500 star cluster. This 
peak in the probability density for the 500 star cluster 
occurs at the highest [^^0]/[^^0] permitted by the 20 
Myr constraint (i.e., the lowest possible mass for SNe II 
progenitors). The other mixing curve is based on the 
lowest [^'^0]/[^^0] ratio consistent with the oxygen iso- 
topic composition of the ISM 4.6 Gyr before present (i.e., 
the present-day [^*0]/[^'^0] ISM, see above). The mix- 
ing curves demonstrate that smaller clusters, represented 
here by our statistical representation of clusters of 500 
stars, were more likely sources of oxygen isotope enrich- 
ment for the solar system than larger clusters (e.g., the 
5,000 star clusters represented by grey contours). 

Results of this analysis show that oxygen ejecta with 
the composition required to explain the oxygen isotope 
ratios of the solar system relative to the ISM 4.6 Gyr 
before present is more than twice as likely to have come 
from a star cluster of several hundred stars than from a 
cluster of several thousand stars. Larger clusters tend to 
produce oxygen too low in [^*0]/[^®0] while considerably 
smaller clusters produce too few SNe II. In our analysis, 
19% of clusters of 500 stars produce oxygen isotope ratios 
suitable to explain the solar composition by supernova 
enrichment. For comparison, the corresponding prob- 
ability for a cluster of 5,000 stars is 7%. Relationships 



betwee n th e mxing curves constrained by the contours in 
Figure 11 the oxygen isotopic composition of the ISM, 
and the oxygen isot opic composition of the solar system 
are shown in Figure [T2| 

The mass fraction of exogenous oxygen in the solar 
parental cloud ca n be estimated from the mixing curves 
in Figures 11 and 12 The masses of ejected o xyg en r ep- 

" lland[l2| 



resented by the two mixing curves in Figures 
are 0.7 and 7.5 Mq for the high and low [ 
cases, respectively. The minimum mass of the enriched 
proto-solar cloud material, Mmc , can be calculated from 
the assumption of a 100% efficiency of injection and si- 
multaneous solution of the mass balance equations 
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(14) 



(15) 



where Cno,Q and Cisq q are the solar concentrations 
by mass of the specified isotopes, Mitq (-^i^cMc) 
is the unknown initial mass of "'^''O (^^O) in the cloud 
material , and A/itq^ej (-^iso.ej) is the calculated mass 
of I'^O (1*^0) added to the cloud by the SNe II ejecta 
constrained by the mixing curves. With the additional 
constraint that the initial [^^0]/[^^0] of the molecular 
cloud was on the Galactic line, such that A/itq^mc — 
Mi8o,Mc(17/18)/4.1, we obtain A^mc ~ 700 to 5,000 
AIq. The mass fraction of SNe II oxygen of 0.1% (e.g., 
100 X 0.7/700), together with SNe II oxyge n mass frac- 



tions 5 to 10 times solar (Rauscher et al. J2002j), suggests 
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a total contribution of SNe II ejecta to the proto-solar 
molecular cloud of ^ 1% by mass. 

4.3. Collateral Isotopic Effects 

Enrichment of the proto-solar cloud by SNe II ejecta 
may help to explain other puzzling aspects of the isotopic 
composition of the solar system. We consider here the 
implications for the isotopes of silicon, carbon, and the 
short-lived radionuclides ^°Fe and ■^^Al as well as for the 
interpretation of pre-solar oxide grains. A caveat to any 
such calculation that attempts to use a single supernova 
source for numerous isotopic and elemental systems is 
the prospect for heterogeneous mixing between the su- 
pernova ejecta and the target cloud material. Nonethe- 
less, in order to assess the collateral implications of 1% 
enrichment of the proto-solar molecular cloud by SNe 
II ejecta for other isotope systems, we derive a general 
mixing equation that illustrates explicitly the relative im- 
portance of elemental abundances and isotope ratios in 
producing mixtures of molecular cloud and SNe ejecta. 
Mass balance for the number of atoms of isotope 1 of el- 
ement n-E-i^, between pre-enrichment molecular cloud 
(MC) material, supernova ejecta (SNe) and the final mix- 
ture (solar, 0) can be written in terms of the total atoms 
for each reservoir j, Nj, and the atomic fractions of nu- 
clide El in reservoir j, {xEi)j, such that 



NmC + A^SNc 

^ (a^B^Mca^MC + (2;£;i)sNoa;sNc 



(16) 



where xsNe = A'sNc/(AfsNG + A^mc), the atomic fraction 
of supernova ejecta in the mixture, and xmc is the atomic 
fraction of original molecular cloud material. Recogniz- 
ing that Xmc = 1 — xsNc and taking the ratio of equation 
( 16 ) for two isotopes of element E we obtain 



2/1 D (^£2)0 



(a^Bi)© 

a;sNc(a;£;2)sNc + (1 ^ a;sNc)(a;_E2)MC 
XSNe{xEi)sNc + (1 - a;sNc) (a;_Ei )mc 



(17) 



Equation (17) is rearranged to obtain an expression for 
the atomic traction of supernova ejecta in terms of el- 
emental abundances and isotope ratios for supernova 
ejecta, pre-enrichment molecular cloud, and the final so- 
lar mixture: 



XSNc 



(18) 



where 



- ^^^Rq ((a;Ei)sNG 



(a;Ei)Mc) 



2/1 



Rmc 



(19) 



In practice we can equate atomic fractions with atoms 
per hydrogen in applying equation (18 1. We will assume 
that the relative abundances of the elements (as opposed 
to isotope ratios) in the precursor cloud were indistin- 
guishable from solar values. Numerical experiments in 
which fictive molecular cloud elemental abundances were 



used confirm that deviations from this simplifying as- 
sumption have negligible effects on the results. 

4.3.1. Carbon isotope ratios 
The solar system [^^C]/[^^C] of 89 is gre ater than the 



typic al present-day local ISM value of ~ 68 ( Milam et al. 
2005|). However, the majority of this difference is ex- 



plained by the increa se in ^^C relative to ^^C by GCE 
over the last 4.6 Gyr ( Prantzos et al.|1996 ). Application 
of equation ([l8| using the sl9 SiNe li model of RHHW02 
([i2C]/[i3C] = 193.9) shows that addition of 1% by mass 
of SNe ejecta from B stars would have raised the pre- 
enrichment molecular cloud [^^C]/[^'^C] from 85 to the 
solar value of 89. There appear to be no observational 
constraints against such a shift. 

4.3.2. Silicon isotope ratios 

A long-standing problem has been an apparent excess 
in ^*Si in the solar system relative to expected values. 
An excess in ^^Si over ^^Si and ^°Si in the solar system 
relative to the ISM 4.6 Gyr before pre sent is evident by 
cornpariso ns with pre-solar SiC grains ( Alexander & Nit- 
tler|1999" ). Mainstream SiC grains come from AGB stars 
that predate the Sun by hundreds of millions to billions of 
years. Galactic chemical evolution should therefore have 
resulted in the younger Sun having greater p^Si]/[^*Si] 
and p°Si]/p^Si] than these earlier- formed AGB s tars, yet 
solar values ar e lower, not higher, by 11 to 12% (Alexan- 
der fc Nittler][l999) . Type II supernovae expel an excess 
of ^*Si relative to the heavier Si isotopes ( Rauscher et al. 



2002 1 and it has been suggested previously that the so 
lar syste m might have been enriche d by ^*Si from su 
pernovae (Alexander fc Nitt ler||1999[ ). Models for SNe II 
ejecta tend to produce too much ■^'^Si relative to ^^Si com- 
pared with the silicon isotopic compositions necessary to 
explain the relationship between sola r sy stem and main- 
stream SiC pre-solar grains (Figure 13), but this may 
be a proble m of spuriou s overproduction of '^"Si in the 
calculations (|Alexander & Nittler 1999). Using equation 



(18), we find that addition of 1% by mass of ejecta from 
the sl9 model of RHHW02 decreases [29Si]/[28Si] and 
P"Si]/[^*Si] in the pre-enrichment solar system molecular 
cloud by 12% and 11%, respectively, making the ISM 4.6 
Gyr before present greater in [^SSiJ/p^Si] and pOSi]/[2SSi] 
than the majority of pre-solar mainstream SiC grains. 
Therefore, the same enrichment process that explains the 
aberrant [^*0]/[^^0] of the solar system could also ex- 
plain most of the excess in ^*Si in the solar system (Fig- 
ure |13[ ). A caveat is that the oxygen isotopic composition 
of the sl9 supernova ejecta model is not ideal for explain- 
ing the anomalous [^^0]/[^^0] of the solar system as it 
is slightly lower in [-'^®0]/[^^0] than progenitors depicted 
with the mixing cuves in Figures [Til and |12l and would 



belt 



therefore require the ISM 4.6 Gyr before present to have 
been higher in [^^0]/[^^0] than today (violating expec- 
tations from GCE). This problem is not severe, however, 
given the uncertain dispersion in ISM [^®0]/[^^0] along 
the Galactic slope-1 line in three-isotope space. 

4.3.3. Short-lived radionuclides 

Assessing the effects of SNe II enrichment on the abun- 
dances of short-lived radionuclides is complicated by the 
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Fig. 13. — Silicon three-isotope ratio plot in which p^Si]/p**Si] 
(ordinate) and [^"Si]/p*Si] (abscissa) are normalized to solar ra- 
tios of 0.05078 and 0.03347, respectively. The plot shows the rela- 
tive positions of the Sun (0), the pre-solar mainstream SiC grains 
(grey circles, Ernst Zinner, pers. comm.), the various supernovae 
ejecta models of RHHW02 (open circles), and supernovae ejecta 
models of NTUKM06 (solid circles). SNe models are labeled with 
progenitor masses (A/q). Results of a mixing calculation yielding 
the composition of the interstellar medium 4.6 Gyr before present 
(black triangle) based on mixing 1% by mass of ejecta from the sl9 
model of RHHW02 with the ISM to produce the solar composition 
is shown with the solid line. The dashed line shows the trajectory 
of the mixing that would be produced using the 13Mq progenitor 
SNe ejecta of NTUKM06. 

opportunity for decay of these radioisotopes in the inter- 
val between SNe ejection and formation of the first soHds 
in the solar system (the so-called "free decay" time). In 
these cases SNe II ejecta isotope ratios ^^^R are mod- 
ified to account for decay of radionuclide i such that 
Vii? ^ */Ji?sNcii exp(-AVr) where T is the mean life of 
nuclide i and At is the free decay time. Gounelle et al. 
([2009;) have argued that a process of enrichment from an 
earlier generation of star formation can explain the short- 
lived radioactive nuclides that are evinced in the solar 
system. Our oxygen isotope results place constraints on 
possible progenitors of these nuclides. For example, ad- 
dition of 1% by mass of SNe II debris from one or more 
BO to Bl stars, as suggested by the [^^0]/[^''0] data and 
represented by model sl9 of RHHW02 ([*^°Fe]/[5*5Fe] = 
8.3xl0-^ pFe]/[H] = 2.61x10-*), would yield the 
amount of ^°Fe (mean life = 2.2 Myr) in the solar sys- 
tem (pOFe]/['5*'Fe] ~ Sx lQ-'^, P^Fe]/fH] = 3.16x10- 5) in- 



ferred from meteorites (Gounelle & Meibom 2008) with 



a free decay time between nucleosynthesis and incorpo- 
ration to solar system precursor cloud material of ~ 11 
Myr. By this time shorter-lived ^^Al and **Ca, t wo nu- 
clides that appear to correlate in the solar system (lSa hij-| 



pal et al.|1998p , would have decayed away, requiring that 
the origin ot these species was separate from that of the 
anomalous abundances of the stable isotopes and ^°Fe. 

4.3.4. Presolar grains 

The oxygen isotope ratios of pre-solar oxide grains 
found in i neteorit e s bea r on the question of the GCE 
of oxygen. Nittler ( 2009 ) argues that models for the ori- 
gins of these grains are consistent with [^^0]/[*''0] equal 



to the solar value of 5.2 for all of the AGB stars from 
which they derive . In all cases, however, fundamen- 
tal problems arise with interpretations of pre-solar grain 
oxygen isotope ratios when the solar system [^®0]/["'^^0] 
is not equal to t he interstellar medium (ISM) today (e.g., 
Alexander & Nittler 1999) . Nittler and colleagues have 
interpreted the pre-solar oxides as having come from a 
range of masses of AGB stars from about 1.2 to 2.2 Mq. 
This corresponds to a sampling of Galactic oxygen over 
5.5 Gyrs based on the lifetim es of the progenitor s stars 
between 6.6 and 1.1 Gyr (e.g.,|Schaller et al.|1992t. That 



is to say, a 1.2 Mq star that enters the AGB phase of 
evolution at the time the solar system was forming sam- 
pled oxygen from the ISM 6.6 Gyr prior to the formation 
of the Sun, or by any measure of stellar lifetimes, very 
early in the evolution of the Galaxy. A 2.2 Mq star en- 
tering the AGB phase at the time the Sun was forming 
represents a sampling of oxygen from 1.1 Gyr prior to 
the formation of the Sun. Therefore, if the interpreta- 
tion of the pre-solar grain data is taken at face value, the 
ISM was characterized by a solar [*^0]/[*^0] of 5.2 for 
7 billion years leading up to the formation of the Sun, 
then sometime between 4.6 Gyr ago and now, the bulk 
[i*0]/["0] of the ISM changed by nearly 30% after hav- 
ing been constant for the prior 7 Gyrs. Such a change 
might be expected once the Milky Way ages to the point 
where high-mass stars are no longer being made, but 
that is not the case yet. Therefore, we conclude that 
it may be necessary to modify models for the origins of 
pre-solar grains. We note that the pre-solar grain oxy- 
gen isotope data cluster about a Galactic [*^0]/[*''0] of 
about 4.1 (Figure 14), suggestive of a causal relation- 
ship between the peak in the pre-solar grain data and 
the Galactic [^*0]/[*^0] (though a coincidence cannot be 
discounted) . In addition, there are pre-solar oxide grains 
that fall below the solar [i80]/[i^0] line (with higher 
than solar [*^0]/[^'^0]) that are not consistent with AGB 
or supernova predictions and so are as yet unexplained by 
any model, suggesting that the origin of high [^*0]/[*^0] 
grains may not be understood in general. 

4.4. Probability of Clusters with Suitable Oxygen 
Isotope Ratios 

The peculiar nature of the solar [*^0]/[^^0] suggests 
that the type of chemical and isotopic enrichment during 
star formation proposed here is not the norm. The like- 
lihood for the enrichment by B stars but not O stars can 
be estimated from the probability for the occurrence of 
clusters of stars having the requisite number of stars (ca. 
500 as opposed to several thousand) and the fraction of 
those clusters that produce supernovae. The cluster mass 
function in terms of the distribution (/) of cluster sizes 
(specified by the number of stars comprising the clus- 
ter, TV*, rather than mass for consistency with our sta- 
tistical analysis) is characterized with a power law (e.g. 
Elmegreen fc Efremov|1997 Parmentier et al.|2008 ) such 
that 



n: 



(20) 



The distribution function is therefore f{N^,) cx Nl~". 
The likelihood of a cluster having from A^min to iVmax 
stars is given by the integral of the distribution function: 
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P 
6 



1.1 Gyr 

(2 2 M,^,^^) 




6.6 Gyr 



Bivariate Probability Density 

(max - 0.D4, contour interval - 005) 

■''0/'^0 = 4.1 

Solar System 



6'°0' = 10^ ln(("0/"'0)/('°0/^''0),s„) 

Fig. 14. — Probability contour plot of the oxygen isotopic compo- 
sitions of pre-solar oxygen-bearing grains (data compilation from 
Larry Nittler, pers. comm.). The peak in the density of pre- 
solar oxygen isotope ratios coincides with the typical Galactic 
["'^^O]/^"'^ O] value of 4.1, represented by the line in the figure. This 
["'^^0]/["'^^0] is significantly lower than the solar value (solar system 
is shown as open circles). Models (Nittler 2009) for the masses 
and ages of the AGB sources of these grains suggest that they 
span many billions of years of Galactic evolution. 



P, 



AN, 



7V„. 



(21) 



N„ 



where a is t he normalizing factor. The value for a in 
equation ( 20 ) that characterizes t he cluster mass function 
is generally round to be nea r 2.0 (Elmegreen & Efremov 
1997| [Haas fc Anders||2010| ). In this case o = 0.109 for 
a total range in TV* from 100 to 10^. H owever, reported 
values for a vary from about 1.8 to 2.4 (Haas & Anders 
20101 with a varying accordingly. For clusters ot lUO to 
1,000 members and a = 2.0, Pan, = 0.25. For a ranging 
from 1.8 to 2.4 the range in Pan, is 0.11 to 0.62. 
The mass generation function used here (equation 



(10)) yields a best-fit relationship between the number 
oFstars comprising a cluster and the fraction of those 
clusters that produce supernovae within the 20 Myr 



timescale (Figure |9|: 



(22) 



-'^SNc = 1 - exp(-6A^*) 

where b = 1.537 x 10"^. From equations ([2l| and ([22]), 
the joint probability for occurrences of SNe ll within 20 
Myr from a cluster of an appropriate range in A^* is 



Pan,,snc 



Ny 



aN^-" (1 - exp(-6A*)) dN^ (23) 



N„ 



Evaluation of equation (23) yields Pan,,snc = 0.10 for 
a = 2.0 with a range due to uncertainty in a from 0.05 
(a = 1.8) to 0.23 (a = 2.4). Our 5% to 20% probability 
is consistent with previous estimates for the fraction of 
stars formed unde r the influence of massive stars adjacent 



molecular clouds (Mizuno et al. 
We must also consider. 



2007 Hennebelle et al. 



2007f. 

cluster in this size range that produces supernovae has 



lowever, that not every 



the right oxygen isotopic composition because ocassion- 
ally some of these clusters do produce the larger O stars 
that are too low in [^®0]/[^^0] to be consistent with en- 
richment of th e proto-solar cloud. The simulation shown 
in Figure [TT] suggests that roughly half of the clusters 
of ca. 500 stars that produce supernovae also have B 
stars as their maximal stellar mass and so produce oxy- 
gen isotope ejecta suitable to explain the solar system 
[^®0]/[^''0]. The probability derived above should there- 
fore be halved, yielding an estimate of 2.5 to 10% for 
the occurrence of a suitable star cluster. These conclu- 
sions are invalid if SNe II oxygen isotope yields are closer 
to the high-[i80]/[i60] values calculated by NTUKM06 
rather than the values reported by RHHW02, LC03 and 
WH07, for example. 

4.5. The Potential for Enrichment by WR Winds 

Another source of oxygen to consider in a star-forming 
region is that produced by the prodigious winds of Wolf- 
Rayet stars (in particular WC stars with He burning 
products exposed at the surface). We conclude that 
this is a less likely alternative to SNe II ejecta. Firstly, 
WR stars evolv e from more massive O stars generally 
( Crowther||2007 ) , making them by their very nature less 
common than less massive B stars that end their lives 
as type II supernovae. Indeed, the ratio of the rates of 
occurence of type II SNe to the rates of occurrence of 
type I b/c supernovae, the lik ely endpoint for WR stars, 
is about 5 (Vanbeveren 2005). Secon dly, the WR phase 
of ev olution lasts for lO'^ years (e.g. Maeder & Meynet 
1994 and references therein) and so the chances ot catch- 
mg winds from this phase of evolution alone, without 
also capturing the ensuing collapse supernova debris, are 
small. Thirdly, rates of mass loss from WR stars are of 
order 10^* to 10~^ Mq /yr, and these rates multiplied 
by the duration of the WR phase of evolution yi elds of 
order one s olar mass of total wind material (e.g. Binns 
et al.|2006[ ). The mass of oxygen released will be consid- 
erably less than the total mass in these winds. Therefore, 
the mass of oxygen released is << 1 Mq compared with 
the minimum of approximately 1 Mq of oxygen liberated 
by low-mass type II supernovae. When considering the 
brevity of the WR phase of evolution, the relative rarity 
of O star WR progenitors, and the relatively low oxygen 
yields from the winds, it seems that enrichment of a star 
forming region in oxygen isotopes from WC stars is not 
as likely as enrichment by SNe II ejecta. 

5. CONCLUSIONS 

We conclude that the comparison between solar and 
extra-solar oxygen isotope ratios in YSOs provides a sig- 
nature of enrichment of the proto-solar molecular cloud 
by supernovae from an earlier generation of star forma- 
tion. Proximity to a previous generation of star forma- 
tion implies that the Sun formed in a molecular cloud 
complex where multiple generations of star formation oc- 
curred. The source of exogenous oxygen that enriched 
the solar system parental cloud is more than twice as 
likely to have come from a smaller cluster of several hun- 
dred stars than from a larger cluster of several thousand 
stars. Tighter constraints on the oxygen isotopic com- 
positions of ejecta from type II SNe would considerably 
improve our ability to use isotope ratios to infer the birth 
environment of the solar system. 
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APPENDIX 



SOLAR SYSTEM [^^0]/[^'^0] 

Specifying a single [-'^^0]/['^^0] ratio for the solar system has three sources of uncertainty. One is systematic errors 
associated with measuring absolute ratios of isotopes. A second is the dispersion in this ratio within the solar system. 
Yet another is the difficulty in obtaining oxygen isotope ratios for the Sun itself. We address these sources of uncertainty 
here. 

Absolute ratio: The [^*0]/[^^0] of standard mean ocean water (SMOW), the commonly accepted reference for oxygen 

(|1976|) to be 0.00205.2 ± 0.0000045 (la), corresponding to a [^^0]/\^^0] 



of 498.7 ± 0.1 
D^^O and H^^o 



isotope ratios, was determined by Baertschi 

This value has remained undisputed and was determined mass spectrometrically using mixtures of 
^^ 0]/[^^0] for SMOW i s less precisely known. Values over the years have generally ranged 



The 



from 0.000365 to 0.00040 ( Santrock et al. 1985), corresponding to [i60]/[i^0] values of 2739 to 2500. The range in 



[i80]/["0] of SMOW based on the well-defined [i^0]/[i80] of SMOW and the full range in reported SMOW [i60]/["0] 
values is from 5.0 to 5.5. Recently, the [i''0]/[i'^0] of SMOW was determined relative to an absolute ["C]/[^^C] for 



the so- called Pee Dee Belemnite standard using measurements of CO2 isotopologue ratios ( Assonov & Brenninkmeijer 
2003b I . The Pee Dee Belemnite (PDB) is a marine fossil composed of calcite that serves as the primary reference 
for terres t rial ca rbon isotope studies. The PDB [I3c]/[i2c] used for this study is 0.011237 ± 0.000030 as determined 
by ICraigl (|1957|) , corresponding to pG]/pC] = 89.0 ± 0.2. The ["0]/[^^0] SMOW value determined by jAssonov fc 
Brenninkrneijerl ( |2003b| ) tied to the |Craig| ( |1957| ) PDB carbon isotope ratio is 0.0003867 ± 0.0000009, or a p^OjTf^ 
of 2586 ± 6. Using this latest value for the absolute abundance of ^^O, the [i80]/["0] of SMOW is 5.18 ± 0.01 (Icr). 
However, the PDB carbon isotope ratio has an associated uncertainty that affects the determination of [^®0]/["'^^0 



for SMOW (Assonov & Brenninkmeijer 2003a|b ). An alternative to relying on a reference car bon isotope ratio is to 



Valkiers et al. (20071 



determine [-'^'OJ/[^''OJ and [^■^CJ/[^^CJ simultaneously for a well-characterized sample of CO2 
performed this experiment on CO2 gas generated by acid digestion of a commonly used carbonate reference material 
NBS-19 calcite. For this CO2 gas these authors obtained [i^0]/[i60] = 0.00038014± (48), [i®0]/[i60] = 0.00208365 ± 
(48), and [^^GJ/f^^C] = 0.0111593 ± (16). Because the difference in carbon isotope ratios between NBS-19 and PDB 
is precisely known, the NBS-19 GO2 carbon isotope ratio corresponds to a [1^G]/[12G] for PDB of 0.0111376 ± (16), 
or [13G]/[12G] = 89.79 ± 0.02 . The differences in oxygen isotope ratios between this NBS-19 GO2 gas and SMOW 



are mass dependent, 
relationship 



It is therefore possible to back out the 0]/[^^0] of SMOW from these new data and the 



NBS - 19, CO, 



^'''-'/ NBS - 19, CO, 



- ITQN 

. if'O/SMOW 



V 



-I'^O/SMOW 



(Al) 



where the exponent (3 characterizes the mass-dependent fractionation law. Values for that relate mass-dependent 



variations in [^^Oj/i^^OJ to those of [^^Oj /j^^O] range from 0.5 1 to 0.53 (Young et al. 2002). For any value for /S 



Al 



the data from iValkiers et al 



(2007), and the accepted [^»6j/[^''0] for SMOW yields a 



2685), well within the range of previous estimates. The resulting 



in this range, equation 
SMOW [i^0]/[i60] value of 0.0003725 {[^^O] 
[i80]/["0] for SMOW is 5.38. 

One can conclude from the discussion above that the absolute [^*0]/[^^0] of SMOW, the primary oxygen reference 
for the solar system, is 5.2 based on the commonly accepted carbon isotope ratio for PDB, and that there is an 
uncertainty in the SMOW [^^0]/[^^0] of about ±0.2 as a result of the uncertainty in the absolute ratio of carbon 
isotopes in the PDB reference material. We show below that the dispersion in the ratio of the rare isotopes of oxygen in 
the solar system is sufficiently small that the uncertainty in the absolute [^^G]/[^'^G] for PDB is currently the dominant 
source of uncertainty in [^'^0]/[^^0] for solar system rocks. 

Dispersion of oxygen isotope ratios in the solar system: Deviations in [^®0]/[^^0] from the value for SMOW arise 
as a result of mass-dependent fractionation of the oxygen isotopes. Most physicochemical processes result in mass- 
dependent partitioning at some level. At issue is how much fractionation exists in solar system materials. Figure |15| 
shows the full range of oxygen isotope ratios for solar system rocks. The most fractionated ratios with the highest S^^O 
are the result of low-temperature water-rock reactions and/or evaporation and are not representative of primordial 
solar system oxygen isotope ratios. Nonetheless, it is clear that all of the data for solar system rocky bodies fall within 
0.1 of [^^0]/[^''0] = 5.2 regardless of level of fractionation. This dispersion in [^^0]/[^''0] is less than the uncertainty 
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Fig. 15. — Three-isotope plot relative to standard mean ocean water (SMOW) showing oxygen isotope ratios of meteorites and terrestrial 
rocks representing Earth, Mars, asteroid 4 Vesta, and the asteroid parent bodie s represented by all classes of meteorites (grey circles). 
The s ou rces of this compilation o f literature data in c lude [Clayton et al.| ([T977t, [Clayton et al.|p 984|, Clayton ct al. (1991), Claytqn[ 
l |1993| , [Clayton fc Mayedaj fl996| l, [Young fc Russell| l|1998| , | Voung et al.| l [Tggg| , [Kobayashi e t al. {'ZOO'd), and Sakamoto ct al. {'ZUOW 
among others. With two exceptions, all data are results ot fliionnation ot riacks, yielding the greatest precision. I'he exceptions arc tlie 
two extreme values obtained on rare materials by secondary ion mass spectrometry presented by Kobayashi et al.| (|2003 ) and Sakamoto[ 
[et al.| l ["2007[| . Uncertainties in each datum are smaller than the symbols used. These data span the lull range ot compositions, including 
highly Iractionated FUN inclusions as summarized by |Krot et al.| l )2010( l. The result for solar wind from the Genesis mission ( |McKeegan[ 
[et al. 2010b) with 2a error bars is shown for comparison. Also shown for comparison are lines of constant [^*0]/[^'^0] based on a value for 
^OW ot 5.2. 



imparted by ambiguity in absolute isotope ratios, as described above. 

Oxygen isotopic composition of the Sun: It is commonly assumed that the Sun has oxygen isotope ratios somewhere 
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along the slope- 1 line shown in Figure |15[ Recently, this assumption was put to the test by mass spectrometri c 



measurements of oxygen isotope ratios or solar wind returned by the Genesis mission (McKeegan et al. 2010a|b ) 



These preliminary results, after correction for mass fractionation in the concentrator where the samples were taken 
(based on "^^Ne/^^Ne fractionation), yield = -99 ± 5%c(2cr) and ^^^O = -79 ± 9%c relative to SMOW. These 

values define an [-'^^0]/[^^0] of 5.07 ± 0.08 {2a), within error of the value for rocks as described above (Figure 15). 
We conclude that the [^^0]/[^^0] of the solar system as a whole, rather than just the rocky bodies, is well defined at 
5. 2 ± 0.2. 

Ayres et al. ( 2006[ ) report oxygen isotope ratios for the Sun based on measurements of the abundances of CO 



isotopologues in the atmospheric layer above the photosphere cooled by rovibrational transitions in CO. The absolute 
[^^0]/[^*0] and [^^0]/[^^0] ratios are sensitive to models for the thermal structure of the middle photosphere (see 
Table 7 in that work), but all of the derived ratios give [^®0]/[^^0] values closer to 4 than 5. The preferred [^^0]/["'^*0] 
and [i60]/[i'^0] ratios of 1700 ± 220 and 440 ± 6 correspond to a [i*0]/[i^O] of 3.9 ± 0.5 (uncertainties are based on 
standard errors about the means). At face value, these results are in conflict with all mass spectrometric measurements 
of solar system materials, including those of the solar wind. We offer no explanation for the disparity between oxygen 
isotope ratios determined for the COmosphere and those for the rest of the solar system. In view of the larger 
uncertainties associated with the former, and the as yet poorly constrained systematic errors associated with solar 
atmosphere thermal models, we do not consider these solar spectroscopic measurements to be remotely as reliable as 
the mass spectrometric measurements. Further study is warranted in all events. 
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